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Abstract: We investigate the behavior of energy momentum tensor correlators in strongly 
coupled large- iVc Yang-Mills theory at nonzero temperature, working within the Improved 
Holographic QCD model. In particular, we determine the spectral functions and corre- 
sponding imaginary time correlators in the bulk and shear channels, and compare the 
results to recent perturbative and lattice calculations where available. For the bulk chan- 
nel imaginary time correlator, for which all three results exist, lattice data is seen to favor 
the holographic prediction over the perturbative one over a wide range of temperatures. 
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1 Introduction 

Alongside with the realization that the quark gluon plasma (QGP) created at RHIC should 
be described as a strongly coupled liquid rather than a gas of weakly interacting quasipar- 
ticles [1], holographic methods have become a standard tool in making qualitative — and 
in a few cases even quantitative — predictions for heavy ion physics [2] . Perhaps the best 
known example of this is the famous conjecture of a lower limit for the shear viscosity to 
entropy ratio, rj/s > l/(47r) [3], which the QGP appears to almost saturate (see also [4-6]). 
In addition, holography has been used to address many complicated dynamical problems 
out of the reach of conventional field theory methods, such as strong coupling thermal- 
ization and particle production (see e.g. [7, 8] and references therein). Within thermal 
equilibrium, it has finally been observed that the behavior of many bulk thermodynamic 
quantities near the critical temperature of the deconfinement transition, measurable with 
lattice methods, can be reproduced to a very good accuracy using holographic models with 
broken supersymmetry and conformal invariance [9-11]. 

On the field theory side, a major obstacle in the quantitative description of the ther- 
malizing plasma is the inapplicability of lattice methods to real time physics — or even to 
the determination of transport coefficients. In the latter case, discussed extensively e.g. in 
[12], some progress has recently been made in combining lattice measurements of Euclidean 
correlators with perturbative results for the corresponding spectral functions [13, 14]. De- 
spite this, important hydrodynamic parameters such as the shear and bulk viscosities are 
still outside the realm of accurate first principles calculations. Recalling that these quanti- 
ties are readily available in the strongly coupled limit of a class of large- A'^c field theories via 
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the AdS/CFT conjecture, it is not surprising that quite some attention has lately turned 
towards a quantitative comparison of lattice, perturbative and gauge/gravity predictions 
for various (mostly Euclidean) correlation functions. Recent lattice studies of energy mo- 
mentum tensor correlators include at least [15-18], while related perturbative work has 
been performed in [19-23] and holographic calculations in [24-30]. Finally, closely related 
studies of sum rules that the associated spectral functions must obey can be found e.g. from 
[31, 32], while an analytic study of the UV limit of different correlators was performed in 
[33]. 

A particularly interesting comparison of lattice, weak coupling and holographic cor- 
relators was reported in [16]. There, it was found that lattice data for various Euclidean 
Green's functions just above the deconfinement temperature of SU(3) Yang-Mills plasma is 
better described by infinitely strongly coupled = 4 Super Yang-Mills (SYM) theory than 
by a leading order perturbative calculation in the original theory. Even though there are 
indications that the inclusion of further perturbative orders acts in the direction of closing 
the gap between the weak coupling and lattice results [20, 21], it is equally worthwhile 
to attempt to improve the description of the system on the strong coupling side. To this 
end, in [30] we addressed the determination of the shear channel spectral function and the 
corresponding imaginary time and coordinate space correlators in the so-called Improved 
Holographic QCD (IHQCD) model, which exhibits a dynamical dilaton field that has the 
effect of breaking conformal invariance and supersymmetry [9, 10]. The study revealed 
important quantitative effects originating from the loss of conformal invariance near Tc, 
and in addition highlighted the importance of performing similar computations in the more 
complicated bulk channel, where conformal theories (such as A/" = 4 SYM) lead to vanishing 
correlation functions. 

In the paper at hand, our aim is to continue and extend the treatment of [30] by 
performing a detailed analysis of the bulk channel of the IHQCD model, concentrating in 
particular on the spectral function at vanishing external three-momentum and the associ- 
ated imaginary time correlator. The latter quantity is of special interest to us due to the 
recent emergence of the corresponding perturbative and lattice results [17, 21], to which 
we can compare our holographic predictions. In addition to this, we will briefiy revisit 
the shear channel, where a Next-to-Leading Order (NLO) result for the spectral function 
has been determined since the appearance of [30], motivating a reanalysis of the IHQCD 
calculation. 

Our paper is organized as follows. In section 2, we review the setup both on the field 
theory and gravity sides, recalling the most important aspects of the IHQCD model. In 
section 3, we next write down the fiuctuation equations, from which both the shear and 
bulk spectral functions are determined, and in addition explain the most important steps 
of the holographic calculation. Section 4 then reviews the existing perturbative results for 
the quantities of our interest, while sections 5-6 contain our holographic results in the two 
channels. A more comprehensive discussion of the results is finally left to section 7, where 
we also draw our conclusions. 

Our notation follows closely that explained in section 1 of [30] . In particular, with the 
exception of the introductory section 2, we will in the following set the AdS radius C = 1. 
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2 Setup 

2.1 Field theory 

We work within pure SU(A'c) Yang-Mills theory at a nonzero temperature T, defined by 
the Euclidean Lagrangian 

5e = l^^dr I d^x \f^,F^, , F;, = d,Al - d.A^ + fe/"'^^^^^ , (2.1) 
with /3 = 1/r. The energy momentum tensor of the theory takes the form 



V(X) = \6,uF^o.pKp - F^aKa = 0,.{X) + J V^(^) > (2-2) 



where we have in the latter stage separated the traceless part 6uu and the anomalous trace 



0ix) ^ T,, = ^-^f;,f;, , (2.3) 



in which /3(g) denotes the beta function of the theory. 

In this paper, we are interested in correlation functions of the shear and bulk operators 
Ti2 and 9 of the above theory, of which the latter we can furthermore replace by the simpler 
quantity Ta, as correlators of Too are known to reduce to contact terms (see e.g. [17]). This 
implies that the retarded correlators we study obtain the forms 

Gf(u;,k = 0) = -i jd^xe''^'9{t){[Tu{t,^),Ti2{0,0)]) (2.4) 

G^ioj,k = 0) = -ij d^xe^"*0(i)([^ri,(t,x), ^r,-,(0,0)]) , (2.5) 

while the corresponding (zero three-momentum) spectral functions read 

Ps,b{io,T) =lmGf^ki^,k = 0). (2.6) 

The relation between these functions and the associated transport coefficients (the shear 
and bulk viscosities) is finally given by 

= lin, PJ<n , (2.7) 

C = hm . (2.8) 

It is also good to recall that in conformal theories, in which the beta function is zero, we 
have T^^ = 0, and thus vanishing bulk correlators and viscosity. 

2.2 Dual gravity system 

In the IHQCD model of [9, 10], the gravitational system involves a gravity+dilaton type 
action 

^ r . r 4 „ 1 

(2.9) 



levrGfi 



R-^id^f + V{<P) 
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while the background metric takes the generic form 



-f{z)df + (ix^ + 



W) 



(2.10) 



where the radial coordinate z is chosen so that the boundary is located at z = 0. The 
functions (j){z), f{z) and b{z) appearing here are determined from the Einstein equations 



W 

b 
X 



AbW^ 



-b^W, 



j{Wf + lbV), 



-XV bW, 



f = ^fbW; 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



in which the dot denotes a derivative with respect to z, and we have defined X{z) = e'^^^\ 
As the notation suggests, this function is found to be dual to the 't Hooft coupling on the 
field theory side, Ac = g'^Nc, while 

m = (2.15) 

b/b 

is related to the field theory beta function. 

In the vicinity of the boundary, the metric function b{z) is required to satisfy 



b{z) 



— > — 

z^O Z 



(2.16) 



where C is the curvature radius of AdS space. For the other functions, the UV limits 
are obtained from the running of the field theory coupling, while the IR behaviors are 
determined by requiring that the model satisfy the confinement criterion of a linear glueball 
spectrum ~ integer [9] (see also appendices A and B of [30]). The dilaton potential we 
use follows the choice of [34] , reading 



ViX) 



12 
£2 



88 4619 2 \/l + ln(l + A) 
^27^ 729 (l + A)2/3 



(2.17) 



It is constructed to reproduce the small A expansion 

V{X) 



12 
£2 



^ 88, 4619,2 



27 



729 



(2.18) 



where the coefficients in front of A and A^ are determined by matching the holographic beta 
function of eq. (2.15) to the perturbative 2-loop result with the identification A = Ac/(87r2). 
The A^ term in eq. (2.17) is finally multiplied by the factor 



V^l + ln(l + A) 
(l + A)2/3 

to ensure that the model satisfies the confinement criterion. 



(2.19) 
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With the gravity system specified, we can now use eqs. (2.11)-(2.14) to solve for those 
bulk field configurations that exhibit a horizon at z = Zh (i.e. satisfy f{zh) = 0) and thus 
correspond to an equilibrium state of the field theory. This enables us to determine a host of 
equilibrium thermodynamical quantities such as the pressure of the system [10, 35], which 
we match to its leading order counterpart in perturbative large- A^c Yang- Mills theory. This 
leads to one further matching relation 

£S 4JV? , , 

which fixes the last unknown parameter of the model. 
3 The calculations 

To determine the correlation functions of the field theory operators T12 and Ta using 
holography, we follow the steps laid out in [25]. First, we introduce perturbations around 
the background metric of eq. (2.10), 

900 = b^f (1 + e^^oo) , 511 = b^l + eHu) , 512 = eb^Hi2 , 955 = j {I + eH^^) , (3.1) 

where e is a power counting parameter. With these definitions, the perturbation Hn 
becomes dual to the operator ^Tu, while H12 corresponds to the shear operator T12. Ex- 
panding the Einstein equations to first order in e, we then find that the metric fiuctuations 
must satisfy the equations 



12 — 


0, 




fx 


p " 


fx 



(3.2) 



iiu + 4- log{b^fX^)Hn + ( ^ _ = , (3.3) 



dz 

where we have set the corresponding three-momentum k to zero and defined 

X{X) = ^. (3.4) 

These equations are to be solved using purely infalling boundary conditions at the horizon, 
most conveniently implemented via an analytic expansion around z = Zh, 

Hs/biz ^ Zh) = {z- Zhr/f'^[l + di{z- Zh) + d2iz- Zh)^ + ...], (3.5) 

and by normalizing the solutions to unity at the boundary. With these initial conditions, 
the equations of motion (3.2)-(3.3) are finally straightforwardly solvable within Mathemat- 
ica. 

With the metric fiuctuations at hand, a standard recipe provides us with simple ex- 
pressions for the shear and bulk spectral functions, 

Ps{oo,T) = .^^JbHmHuHl^ , (3.6) 
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where the combinations on the right hand sides can be shown to be independent of z. 
In practice, their evaluation however turns out to be most convenient near the boundary, 
leading us to the expressions 

«(..T) = 6X|5m_Ji_, (3,9) 

where = and s(T) = h\l{^G^ stands for the entropy of the model. These are 

the quantities we will be inspecting in the following sections. 

4 Perturbative limit 

Before proceeding to the results of our holographic calculations, let us briefly review what 
is known about the behavior of the shear and bulk spectral functions in weakly coupled 
SU (iVc) Yang-Mills theory. This is helpful in particular for the analysis of the UV (large-w) 
behavior of our results, as due to asymptotic freedom all physical correlators are expected 
to reduce to their perturbative limits as a; — )■ cxd. 

At the moment, the k = spectral functions of both channels are known up to and 
including their respective NLO terms in perturbation theory [21, 23]. In the shear case, we 
can read off the result from eq. (4.1) of [23], obtaining (note an additional factor of -1/16 
due to differing definitions of the shear operator) 

= — (l + 2n.)|l-^f^ + c/>^(^^^)')|+0(A^) (4.1) 

^ u:\ (4.2) 



aj->oo IGOvr 

where (Ia = N"^ — 1, = l/(e^/"^ — 1), and (j)j^{uj/T) is a numerically evaluatable dimen- 
sionless function that behaves like T^/lj^ in the — )• 00 limit. It should be noted that this 
result misses a number of terms proportional to lv 5{u}), which give important contributions 
to the shear sum rule but are only known to leading order, cf. e.g. eq. (4) of [18] . 

In the bulk channel, the perturbative spectral function consistent with our earlier 
definitions is obtainable from eq. (4.1) of [21]. Multiplying this result by 1/9 and choosing 
the constant cq as g'^CQ = ^^^"^ , where /3(Ac) is the beta function of Yang-Mills theory, we 
obtain 

dA 576TT \l ^ 2 M 87r2 V 3 3 




+ 0{\t) (4.3) 

(4.4) 



324(47r 



where 4'^{uj/T) is again a numerical function, whose behavior was analyzed in quite some 
detail in [21]. An important difference to the shear channel result is clearly the appearance 
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Figure 1. Left: the IHQCD shear channel spectral function displayed in the region of small 
frequencies for three different temperatures. The dashed curve represents the large-w limit of the 
SYM result, cf. [29]. Right: the behavior of the spectral function at large frequencies for the case 
of r = STc- The black curve stands for our IHQCD result, while the two dashed red lines denote 
the NLO perturbative result evaluated with two different renormalization scales [23]. 

of the 't Hooft coupling in the leading large-w behavior of eq. (4.4). Together with the 
realization that the renormalization scale, with which the coupling runs, is in the limit 
u ^ T necessarily proportional to w, this implies that the leading UV behavior of the bulk 
spectral function takes the form of a T-independent constant times cj^/(ln w/A^jg)^. In 
addition, one should note that in the bulk channel, no terms of the type u}5{uj) appear at 
least at the orders considered above. 

When analyzing the perturbative results for both the shear and bulk spectral functions, 
an important thing to note is that even at high temperatures — and thus weak coupling 
— the above expressions are not valid in the limit of very small uj. This is due to the 
multitude of soft scales that enter the calculation at small momentum exchange and require 
complicated resummations to be performed when entering the regions of oj of order gT, 
g^T and ultimately g'^T (see e.g. [37]). While with the Hard Thermal Loop resummation 
performed in [21] the above bulk result should be correct down to momenta of order oj ~ gT, 
it is clear from the plots of [21, 23] that in both the shear and bulk cases, the perturbative 
results begin to lose accuracy when uj <T. In particular, this implies that even the leading 
order transport coefficients are not available from the above expressions, and that when 
comparing our holographic results to them, one should only expect quantitative agreement 
at w » r. 

5 Holographic results in the shear channel 

The shear spectral function was first determined within IHQCD in [30], and is reproduced 
for vanishing external three-momentum in fig. 1. Comparing to the corresponding result 
in the conformal A/" = 4 SYM theory [24, 29], whose asymptotic behavior is represented by 
the dashed blue curve on the left, we see that the effects of conformal invariance breaking 
are largest near the IHQCD deconfinement temperature Tc, while already at T = lOTc the 
IHQCD result is rather close to the SYM one. In the — t- limit, each of the curves 
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Figure 2. The bulk spectral function evaluated for a set of different temperatures in the region of 
small Lo. The temperature dependence of the bulk viscosity can be read off from intercepts of the 
curves at cj = 0. 



furthermore reproduces the well known result of rj/s = l/(47r), equally valid in IHQCD as 
in the SYM theory. 

Proceeding to larger values of w, we display the behavior of the shear spectral function 
on a log-log scale in fig. 1 (right), where it is further compared with the perturbative result 
of eq. (4.1). This reveals a clear discrepancy between the two results in the UV region, 
which one can understand using the analytic WKB calculation of [30]. With the help of 
eq. (2.20), one namely easily obtains as the limiting behavior of the IHQCD shear spectral 
function 

p.(^,T) ^Lo\ (5.1) 

which deviates from the perturbative limit of eq. (4.2) by a factor of 4/9. This should not 
come as a surprise considering that IHQCD is only a two-derivative model, but nevertheless 
highlights its limitations in describing the UV dynamics of the physical theory. 



6 Holographic results in the bulk channel 

The bulk channel has been studied in various conformality breaking holographic models 
already in several works [25, 26, 36], but typically concentrating only on the behavior of 
the bulk viscosity to entropy ratio as a function of T. In this section, our aim is to extend 
this treatment to the evaluation of the full IHQCD bulk spectral function. 

Beginning from the limit of small frequencies, we first display in fig. 2 the behavior 
of pb{u},T) for temperatures ranging from Tc to lOTc (cf. also fig. 1 (left)). In accordance 
with our expectations, we observe a decrease in the values of the function with increasing 
T, signifying the approach of the system towards the conformal limit. From the intercepts 
of the curves at w = 0, one can furthermore read off the values of the bulk viscosity C at 
different temperatures, leading to a behavior consistent with that shown in fig. 8 of [26]. 
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Figure 3. Left: the behavior of the IHQCD bulk spectral function shown over a wide range of 
frequencies at three different temperatures. Note that unlike in most of our other plots, uj has here 
been scaled dimensionless by and not T. Right: a comparison of the IHQCD (solid black curve) 
and perturbative (red dashed lines) results for the bulk spectral function. The perturbative result 
is taken from [21]. 



Next, we study the large-ci; behavior of the bulk spectral function in fig. 3. On the 
left, we demonstrate, how the combination pf){uj,T)/uj'^ undergoes a sharp transition from 
a T-dependent behavior at small frequencies towards a T-independent l/(lna;/Tc)^ 

limit at UJ ^ T. On the right, we on the other hand specialize to the case of T = 3Tc, 
displaying the holographic result together with the perturbative one, eq. (4.3). This fig- 
ure demonstrates a remarkable fact: not only is the form of the asymptotic l/(lnw/A^|g)-^ 
behavior (note that Tc ~ ^ms) of perturbative result reproduced by our IHQCD cal- 
culation, but even the overall coefficient in eq. (4.4) appears to agree with our numerics. 
We find this possibly coincidential fact very surprising, considering the missing 4/9 factor 
encountered in the shear spectral function. We see no a priori reason, why corrections from 
higher derivative terms in the holographic action should be present in the shear channel 
but absent from the bulk one. 

Ideally, it would of course be pleasing to be able to derive the logarithmic UV behavior 
of the bulk spectral function analytically, following a WKB expansion similar to that 
performed in the shear channel in [30]. In the bulk case, this, however, turns out to be 
a rather demanding task due to the appearance of logarithmic terms in the fluctuation 
equation, originating from the z — )• limit of the quantity 

3X = ^ = ^^J-. (6.1) 
A dlogb z->o logz 

Scaling the radial variable according to z — )■ z' = ujz, the large-a; (i.e. small-z) limit of 
eq. (3.3) namely becomes (with A denoting an arbitrary scale parameter) 

ffii + l --fl+ Q., A/^2 V li \ / , \hii + Hu = 0, (6.2) 
I z\ 9(log Az/o;)^/ z| log Az/wl I 

which without the logarithmic terms would lead to the usual behavior of the spectral 
function. In the presence of the conformality breaking w-dependence, solving the equation 
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Figure 4. The bulk channel imaginary time correlator of eq. (6.3) evaluated for two different 
temperatures in IHQCD (solid black curves) and perturbation theory (red dashed curves), and 
compared with the lattice data points of [f7]. The perturbative result is again taken from [21]. 

however becomes much harder, and in particular leads to the numerically verified appear- 
ance of logarithmic suppression in the spectral function. It is in any case worth noting that 
the resulting l/(ln w/Tc)^ behavior of pi,{u},T) enters eq. (3.9) solely through the z — )• 
limit of Hii{z), and not via the factor X'^ ~ /3(A(z/i))^/A(2;/i)^. This is in clear contrast 
with the perturbative limit in eq. (4.4), in which the logarithmic behavior is due to the 
running of the gauge coupling. 

Having the bulk spectral function now at hand, a natural application is clearly the 
determination of the corresponding imaginary time correlator, for which both perturbative 
and lattice results exist. To this end, we plug our function pi,{uj,T) to the relation 

G{t,T)= —pb{u:,T) { ' , P^l/T, (6.3) 

Jo 

obtaining the result displayed in fig. 4. Our holographic prediction is seen to agree with 
the lattice data better than the weak coupling result over a wide range of temperatures, 
the difference being (not surprisingly) most pronounced close to Tc. 

Finally, a different way of inspecting the imaginary time correlator is to look at its 
value at the symmetry point r = 1/(2T) as a function of temperature. This we do in fig. 5, 
where G{t = 1/{2T),T) is displayed, normalized dimensionless by T^. The plot indicates 
a rapid decrease in the quantity as the temperature is raised above Tc, to be contrasted 
with the slow increase of the corresponding quantity in the shear channel, shown in fig. 6 
of [30]. This fact can clearly be attributed to the system approaching conformality in the 
limit of high temperatures. 

7 Conclusions 

In the paper at hand, we have studied finite temperature correlation functions of the en- 
ergy momentum tensor of large- A'^c Yang-Mills theory, concentrating on the shear {T12T12) 
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Figure 5. The imaginary time correlator of eq. (6.3), normalized by and plotted as a function 
of temperature at the symmetry point t = 1/(27"). 

and bulk {TaTjj) channels at vanishing external three-momentum. In particular, after 
determining the spectral functions and associated imaginary time correlators in the Im- 
proved Holographic QCD (IHQCD) model [9, 10], we performed a detailed comparison 
of our results with state-of-the-art perturbative and lattice works. Clearly, the domains 
of validity of the three methods do not always overlap. Perturbation theory requires the 
gauge coupling to be small, which is formally only realized at asymptotically large T or 
uj, while holographic methods work best in the strongly coupled, yet conformal limit, to 
which systematic corrections are accounted for in the IHQCD model. Finally, while being 
a fundamentally nonperturbative first principles method, lattice QCD is unfortunately re- 
stricted to the Euclidean formulation of the theory, and thus only provides results for a 
limited set of observables. 

Comparing the IHQCD shear and bulk spectral functions with their perturbative coun- 
terparts, cf. figs. 1-3, we witnessed an expected pattern, in which conformal invariance 
breaking effects were seen to be largest near Tc, but rapidly decrease with increasing 
temperature. Furthermore, we saw that in the large-w limit of both the shear and bulk 
channels, the parametric dependence of the perturbative spectral functions on u {uj^ and 
a;^/(ln tj/A^jg)^, respectively) was correctly reproduced by IHQCD. A closer inspection 
further revealed that while in the shear channel the a; — )• oo limit of the perturbative re- 
sult was larger than the IHQCD one by a factor 9/4, surprisinly in the bulk channel the 
asymptotic limits perfectly coincide. We find this quite remarkable, considering that the 
nonzero value of the bulk correlator is entirely due to the conformal invariance breaking 
buih into IHQCD. 

For a set of Euclidean quantities — the imaginary time correlation functions — we 
were able to perform comparisons between IHQCD, perturbation theory and lattice Monte 
Carlo results. A direct comparison of the bulk channel correlator Gf,(r, T) was performed 
at two temperatures, 1.65Tc and 3.2Tc. The results showed the lattice data consistently 
prefering the holographic prediction, though at higher temperatures the difference was seen 




T/T, 
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to somewhat diminish. 

Finahy, we note that in a recent paper [38], an IHQCD calculation closely related to 
ours was performed for the correlators of the pseudoscalar operator TrF^i/F^i,. It is in- 
teresting to compare the results reported in section 4 of this paper to ours, concerning in 
particular the asymptotic large-w behavior of the k = spectral function. While pertur- 
bative arguments suggest that the pseudoscalar spectral function should behave similarly 
to our bulk result (cf. ref. [21]), the authors of [38] argue that their numerical data at large 
frequencies is consistent with a pure behavior. This being the case, it would clearly be 
crucial to understand the physical origin of the differing behavior, perhaps by performing a 
WKB type expansion in the asymptotic region of both channels. This calculation, as well 
as an analysis of the correlation functions of the Tr F^yF^y operator, we however leave for 
the future. 
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